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A system of N Brownian particles suspended in a nonuniform heat bath is treated as a ther- 
modynamic system whith internal degrees of freedom, in this case their velocities and coordinates. 
Applying the scheme of non-equilibrium thermodynamics, one then easily obtains the Fokker-Planck 
equation for simultaneous Brownian motion of N particles in a temperature gradient. This equation 
accounts for couplings in the motion as a result of hydrodynamic interactions between particles. 



I. INTRODUCTION 

In a former paper Q, we reconsidered within the context of nonequilibrium thermodynamics for systems with 
internal degrees of freedom Q , the problem of Brownian motion in a temperature gradient and obtained in a simple 
way, for this case, the Fokker-Planck equation for one Brownian particle. This equation had been originally derived 
in 1968 by Zubarev and Bashkirov H using the more elaborate statistical mechanical methods developed by Zubarev 
for treating nonequilibrium systems. 

More recently Shea and Oppenheim [Q applied statistical mechanical methods based on the use of projection 
operators to the study of Brownian motion of one particle in a nonequilibrium bath. They obtain a Fokker-Planck 
equation which, in the presence of a temperature gradient, agrees with the one derived by Zubarev and Bashkirov ||, 
and thus also with the one obtained from the nonequilibrium thermodynamic scheme for systems with internal degrees 
of freedom. The concept of internal degrees of freedom had been introduced already in 1953 into the formalism of 
nonequilibrium thermodynamics by Prigogine and Mazur who pointed out that in this way Fokker-Planck- like 
kinetic equations could be derived. 

In this paper we generalize our previous result to the motion of many spherical Brownian particles suspended in 
a heat bath in the presence of a temperature gradient. In §2 the system of N Brownian particles suspended in a 
nonuniform heat bath is considered as a system with internal degrees of freedom, consisting of the positions and the 
velocities of those particles. As a consequence, the probability density function in the phase space of the N Brownian 
particles plays the role of a thermodynamic density variable of a many component system, for which we give the Gibbs 
equation, containing the Massieu function conjugated to that density variable. Using Gibbs entropy postulate we are 
able to express this Massieu function in terms of the N Brownian particle distribution function. In §3 we discuss the 
form of this distribution function in the state of local equilibrium, in which each particle has a velocity distribution 
in equilibrium with the temperature at its position. 

We are then able not only to state in §4 the conservation laws obeyed by the energy density and the densities 
of the various components with different internal degrees of freedom, but also to calculate the entropy production 
accompanying non uniformities, in ordinary coordinate space, as well as in the phase space of the internal degrees of 
freedom. 

In §5 we follow the principles of nonequilibrium thermodynamics to formulate the linear phenomenological laws for 
the dissipative fluxes occurring in the entropy source strength. These phenomenological laws contain couplings between 
fluxes pertaining to different particles, which if realized find their origin in hydrodynamic interactions. Substitution 
of the phenomenological laws into the conservation equations gives rise on the one hand to the desired Fokker-Planck 
equation for N particles, and on the other to a coupled differential equation for the temperature distribution of the 
heat bath. 

In §6 we discuss the result obtained, also in comparison to - and in the light of - the merits of a more fundamental 
statistical mechanical approach. Moreover we address the problem why and when the inclusion of internal degrees of 
freedom in the formalism of nonequilibrium thermodynamics gives rise to valid and acceptable results. 
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II. THE BROWNIAN GAS AS A THERMODYNAMIC SYSTEM WITH INTERNAL DEGREES OF 

FREEDOM 



Consider a gas of N Brownian particles of mass m suspended in a heat bath at rest with constant mass density pn 
and fixed total volume V. A nonuniform temperature field T(r) exists in the system (Brownian gas and heat bath) 
at position r. 

The Brownian particles are assumed to be essentially hard spheres and to have no further direct interactions. Their 
center of mass positions in the heat bath are Rj, i = l,2,..,N, their velocities U^. The probability density for the 
Brownian particles to be at time iina state defined by a point T = (R N ,XJ N ) = (Ri,R 2 , .., Rat; Ui, U2, .., Us) in 
the 6N dimensional velocity-coordinate space will be denoted by P N (T, t) and satisfies the normalization condition 

3 N 



p"(T,t)dr = 1. (i) 

The quantity NmP N then represents the mass density of Brownian particles having velocities U w and positions TL N . 
The system itself, Brownian gas and heat bath, has mass density 

p(r,t) = P h + p B (r,t) = PH + J dTj2mS(K^r)P N (T,t) = p H + ^(m^R* - r)>. (2) 

i i 

Here /9e(r,i) is the mass density of the Brownian gas; the bracket (..} denotes an average over the distrubution 
function P N (T,t). 

The coordinates Ri and velocities may be considered as 6N intermal degrees of freedom of the system, in the 
sense of nonequilibrium thermodynamics Taking this point of view, which implies that the system can be looked 
at as a system of many components with mass densities NmP N (T, t), the Gibbs equation for the total entropy S of 
the system may be written as 



5S 



J dr^-Spe - Nm J dTQ(T,t)SP N (T,t), (3) 



where e(r, t) is the energy per unit of mass and the Massieu function conjugate to the density NmP N . In writing 
down equation (||) account has been taken of the fact that pjj is constant. 

Before carrying out the program of nonequilibrium thermodynamics (which is to calculate in the first place the 
entropy production occurring in the system as a consequence of the nonequilibrium distribution of temperature and 
internal degrees of freedom) and thus to derive the Fokker-Planck equation for N Brownian particles in the presence 
of a temperature gradient, we need to establish an expression for the Massieu function 0(r), relating it to the 
(thermodynamic) variables P N (r, t). For this purpose we make use of Gibbs entropy postulate and write the entropy 
S whith respect to its local equilibrium value in terms of a distribution over possible internal states 

r P N (r t) 

S = kJ dTP N {T.t)\n 1 ^-l + S l -<: (4) 



Here P l (T) and S Le ' are the distribution function and the entropy at local equilibrium; k is Boltzmann's constant. 
By local equilibrium we mean that each Brownian particle has at its position a distribution over internal (velocity ) 
states in equilibrium with respect to the temperature at that point. The entropy S Le ' , which depends only on the 
temperature distribution and the density distribution of Brownian particles in the heat bath, obeys the standard 
Gibbs equation 

SS = J dr(±5pe - ffs PB ), (5) 

where ps (r) is the chemical potential of the Brownian gas per unit of mass in internal (local) equilibrium at temper- 
ature T(r) and density ps- 

On the other hand we have from expression (||) for the differential of S 



Comparing Eqs. (|||) and (g) we find using also Eqs. (||) and (||) 
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A/m"'" p^.(r) 

with 

Thus in local equilibrium the Massieu function 0' e (T) is uniform in internal velocity space. 

III. THE N PARTICLE DISTRIBUTION FUNCTION AT LOCAL EQUILIBRIUM 

It remains to determine the local equilibrium distribution function P N (T,t). We first note that in full equilibrium 
(in the absence of a temperature gradient) P N is given by the canonical ensemble distribution function 



P e N q = ±exp{(F N -n\!,T\. (<M 



with T the equilibrium temperature of the heat bath, H N the energy of the N Brownian particles and F N the N 
particle free energy 

H N = Y / \mUl (10) 



F N = F N {T,V) = Nmf B (T,V/Nm). (11) 

Here fs is the free energy of the Brownian gas per unit of mass which (for N sufficiently large) is a function of 
temperature and density only. 

In Eq.(|l^) the direct interactions between Brownian particles (even their hard sphere interactions) are neglected, 
as the system is assumed to be sufficiently dilute. 

We shall now take P N at local equilibrium to be of the form 

PL = ■^exp{(mf l -H l )/kT(R l )}, (12) 

with 

Miji 

fi = / s (T(R 4 ), Pe ), Pe = — , (13) 

Hi = \mUl (14) 

where T(r) is a given non uniform temperature field and p e the uniform equilibrium particle density. 

The free energy which is determined by inserting (12) into (1) and using Stirling's approximation, 
In AH = NhiN — N, has the classical form 

fi = In In 1 = 

m m 2 m 

= n B (Ri) - psiR^p- 1 . (15) 

Here ps (r) is the explicit form of the previously introduced chemical potential per unit of mass of the ideal Brownian 
gas and Pb{ y ) — kT(r)p e /m its - perfect gas - pressure. 
Thus P^ e can be written in the form 

^ =11^ — kfW) — } - (16) 
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If expression ( [Tq ) is inserted into Eq.(Q) we obtain for the Massieu function Q(T,t) the alternative form 

e ( r, t ) = ^i„P»(r, t ) + i£^ + ii» W . ,17) 

We are now in a position to carry out the nonequilibrium thermodynamics analysis which leads to the derivation of 
the N particle Fokker-Planck equation. In the next section we shall establish the basic equations necessary for this 
derivation. 



IV. CONSERVATION LAWS AND ENTROPY PRODUCTION. 

Let us calculate the rate of change of entropy by differentiating Eq.([|) with respect to time 

dS f , 1 dpe Ar f mrt/ _ ,dP N (T,t) 

— = / dr— - Nm / dr&(T,t) ^- Lj -. (18) 

dt J T(r) dt J y ' dt v ' 

To proceed with the analysis we need the conservation laws obeyed by the energy density pe and the distribution 
function P N (T,t). 

We shall first consider the latter: in V space P N obeys the continuity equation 

i i 

where the Ju; are components of fluxes in 3N - dimensional velocity space, which account for the interaction of the 
Brownian particles with the heat bath. We have neglected all direct interactions between these particles also in writing 
down Eq.fllSp. This is a reasonable approximation since in our dilute system the particle dynamics are overwhelm- 
ingly determined by their dissipative interactions with the heat bath. It is indeed our aim to use nonequilibrium 
thermodynamics to find the laws obeyed by the fluxes Jui and thus to derive the equation describing the resulting N 
particle random motion problem (also when the heat bath is not uniform). 

Next to Eq.@ we need the law of conservation of energy which can be formulated as follows 



dpe 

~di 



^ = -V-J„ (20) 



where 3 q represents a heat flux defined in the reference frame in which the heat bath is at rest. 

Substituting in Eq. (|lj) the conservation laws (|l9| ) and ( ^0| ) we can write the rate of change of entropy as the sum 
of two terms 



dS d e S diS 
dt dt dt 

The first term is the entropy supplied to the system, and is given by 

d P S f ,„ J 



(21) 



dt=-J dSl -^ ^ 



while the second term represents the total entropy production and is equal to 



-£ = I drY^{mp^(R i )P N (T)J q<i (R i ).VT- l (R i ) - J^-k—ln^}. (23) 



In Eq.(|22|) ft is the surface of the heat bath and dfl a vector of magnitude dH, normal to the surface pointing 
outward. The vectors J 9 ,i in Eq.(p3|) represent modified heat fluxes at the positions of the Brownian particles and 
are defined as 

J gi< (Ri,Ui) = J,(Ri) - MR^UlU;. (24) 

They contain contributions from the kinetic energy flows of the Brownian particles themselves. 
In deriving Eq. (p2|) and (23) use has been made of: 
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1- Gauss' theorem in coordinate space and the fact that the Brownian particles are confined to the heat bath, 



2- the fact that upon partial integration in velocity space the fluxes Ju, vanish as the velocities ± Uj tend to 
infinity; and 

3- the form(s) (@) and/or @ of the Massieu function 6(T). 

The entropy production (23), which is a positive quantity in accordance with the second law of thermodynamics, 
is written as an integral over the N-Brownian particle phase space. We shall now demand, in agreement with the 
extension of nonequilibrium thermodynamics to systems with internal degrees of freedom |^| , that not only the integral 
( p3| ) be positive definite, but also its integrand itself, 

f = I < 25 » 

d P N 

a r = -J2mps 1 (^i)P N ^,i • Z^VIi -^Ju,- k— ln(-^) > 0. (26) 

i i l- e - 

We have introduced here the notation Ti = T(R^). The quantity or, the entropy source strength in T space, is a sum 
of products of dissipative fluxes 3 q i and 3u i respectively, and their conjugate thermodynamic forces, the gradients 
of temperature in coordinate space and of the Massieu function O(T), cf. Eq.(Q), in velocity space. 

V. PHENOMENOLOGICAL RELATIONS AND THE FOKKER-PLANCK EQUATION FOR 

N-BROWNIAN PARTICLES. 

Following the scheme of nonequilibrium thermodynamics Q we shall now establish linear phenomenological relations 
between the fluxes and thermodynamic forces occurring in expression (^6|) for T . Since the system is isotropic these 
relations are 

F) P N 

mp B 1 P N J q , l = -Y^Lt^T^VT, fc£i TiUj — lu(_) (27) 

3 3 3 Le - 

d P N 

Ju, = -^W^rVZ} - k^Ljj^ — ln^) (28) 

3 j 3 Le - 

Assuming Onsager's principle of microscopic reversibility to hold for the system studied^, the coefficients in Eqs. 
(^7]) and (|2^) satisfy the Onsager Casimir reciprocal relations. 

LTiTj = LTjTi, (29) 

LiJiTj — —LtjUi, (30) 



(31) 



1 Onsager's theorem is based on the property of microscopic reversibility which holds for fluctuations, in particular, around 
equilibrium. But in the system studied we would have to consider correlations of particle velocity fluctuations at different 
positions and temperatures around a quasi-stationary nonequilibrium state, for which the theorem does not necessarily hold. 
Neverthless it seems legitimate to use the theorem when one takes into account that the velocity correlations at two different 
position i and j, which find their origin in hydrodynamic interactions, are sufficiently short ranged in space for temperatures 
to be taken equal as in equilibrium. 
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Defining the heat conductivity coefficients Ay , the thermal acceleration coefficients 7^ and the friction coefficients f3, 



1 :i 



L T% T 3 PB 

AlJ ~ mP N Tf [6 ' 



_ L UiTj 



N = -p^' ( 34 ) 



relations (|27|) and (|28|) become taking into account also (J30|) , as well as the fact that in ( p3| ) Tj ~ in good 
approximation (see also note previous page) 

mp B 1 P N 3 q , i =-J2 m\ l3 p B 1 P N VT 1 + £ m lji {P N V j + (35) 
3 j 3 

urn. f)p N 

Ju, = -E rn^P^Tj/Tj - £ m Ai(P W U j + ^^)- (36) 
j j J 

Eqs. ( |35| ) and (^) (as did Eqs.(^) and (p8|)) give expression to the fact that the dissipative fluxes at different (par- 
ticle) positions can be coupled to each other: the coefficients /3y represent the well-known mutual friction coefficients 
of spherical bodies which result (and are calculated) from their hydrodynamic interactions as they move through the 
fluid they are suspended in. Likewise the coefficients 7y and Ay-, i ^ j could arise from hydrodynamic interactions. 
We do not know whether such effects have previously been found. As far as the heat conductivity matrix is concerned 
its dominant contribution to \u is the heat conductivity of the heat bath itself, while all other contributions and 
matrix elements arise from particle motions. 

We shall now first formulate, using Eq. (|35|) , the phenomenological law obeyed by the total heat current in the heat 
bath. For this purpose we multiply both members of ( |35| ) by <5(Ri — r), sum over all particles i, and integrate over 
T-space. We then obtain the following law 

J' 9 (r) = - J A(r,r')VT(r')rfr' + J T (r, v')p B {v')Y B {v')dv' . (37) 



Here the total heat flux J' q at r is given by (c/. Eq.(24)) 



J' 9 (r) = J,(r) - E<^U?UiJ(Ri - r)), (38) 

i 

while Vefr'), is the mean particle velocity, 

V fl (r) = ^(JmWfRi-r)). (39) 

i 

The heat conductivity kernels A(r,r') and 7(1% r') are defined by the following averages 

A(r,r') = mpzHY, M(R< ~ r W R ; " O), (40) 



7 (r,r> S (r')V B (r') = (^ 7ljra U l( 5(R, - r)*(Rj - r')). (41) 

ij 

These kernels, which for r ^ r' result as stated above from hydrodynamic interactions, can be assumed to be 
sufficiently short ranged, so that from a truly macroscopic point of view Eq. (^) has the approximate form 



G 



J' q = — AVT + 1Pb Vb, (42) 



with the macroscopic heat conductivity 



and 



A = f A(r,r')dr', (43) 

7 = J 7 (r,r')dr'. (44) 

The last quantity is related to the thermal diffusion coefficient of the Brownian gas (for the exact relation of 7 to the 
thermal diffusion coefficient see eg. ref. 0). 

Substitution of Eq.([l2|) into the law of conservation of energy ( pfj| ) gives rise to a differential equation for the 
temperature which is coupled to the first few Brownian particle velocity moments (cf. also Eq. d3S|). These moments 
can be calculated from the N particle distribution function P N (T,t), which in turn is the solution of the differential 
equation obtained by substitution of Eq. (pfl) into Eq. (|l9j) 



l ij j j j 

This is the N particle Fokker Planck equation we set out to derive from the scheme of nonequilibrium thermody- 
namics as applied to systems with internal degrees of freedom. In the presence of a temperature gradient a coupling 
occurs, as could be expected, to the value of that gradient at the position of each particle. 

In our derivation of Eq. (|45|), direct interactions between the Brownian particles were completely neglected. It 
must be pointed out, however, that this derivation can easily be extended to the case that direct interactions ipij — 
¥>(|Ri — Rj|) play a significant role. 

To perform that extension the local equilibrium disatribution function ( |l2| ) must be adequately modified replacing 
Eq. © by 

* = \«iUf + ^tpy , (46) 
3 

thus attributing one half of the interaction energy of particle i and particle j to particle i and the other half to particle 
j. It is then not possible anymore to determine fi explicitely in a simple manner. However with fi = [is —"PbPZ 1 an( i 
P^ e therefore given by 



i 3 



(47) 



(which reduces, whithout direct interactions and the consequent perfect pressure value for ps, to Eq. (16)) the analysis 
can proceed as before. One must then take into account the continuity equation obeyed by P N in the presence of 
direct interactions and obtains in a manner completely analogous to the one described above (thereby introducing 
also a new modified heat current (cf. Eq. (p4])) the full Fokker-Plank equation, 

^ + yu..^T m -^—-y-. \0-(P N U ■ + ^ — ) + 7" P N -^} (48) 

i %3 

The Fokker-Plank equation (Q) now remains applicable to systems which are not necessarily dilute. 



VI. CONCLUSION AND DISCUSSION 



In the preceeding sections we considered an ensemble of Brownian particles suspended in a nonuniform heat bath 
as a system with internal degrees of freedom and applied to it the scheme of nonequilibrium thermodynamics. We 
showed that it is then a simple matter to obtain the Fokker Planck equation for simultaneous Brownian motion of 
N particles in the presence of a temperature gradient. In fact the derivation along these lines turns out to be no 
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more involved than the corresponding one for a single Brownian particle, even though many particle hydrodynamic 
interactions have to be taken into account. It is the simplicity of the nonequilibrium thermodynamic approach which 
forms its main attraction, in particular when it concerns, as it does in this case, phenomena which occur somewhere 
on the interface between macroscopic and microscopic physics. 

This thermodynamic approach should be compared to statistical mechanical methods to derive the basic equation 
describing many particle Brownian motion. In an equilibrium heat bath statistical theories where given by Deutch 
and Oppenheim Q and Murphy and Aguirre Q. While the statistical mechanical treatments are more satisfactory 
from a fundamental point of view, they do not seem in the present stage to provide substantially more detailed results. 

It should be noticed that in the application given above the extension of nonequilibrium thermodynamics to systems 
with internal degrees of freedom already contains elements of statistical physics beyond those present when the concept 
was introduced ||. Thus the theory developed here makes use of the probability density function in the phase space 
of the N Brownian particles and of Gibbs entropy postulate to relate the entropy of the system to this distribution 
function. 

In their original paper Prigogine and Mazur stated that it was their object, in considering irreversible processes 
connected to internal degrees of freedom, to extend the domain of applicability of nonequilibrium thermodynamics 
(d'etendre le domaine d'application de la thermodynamique des phenomenes irreversibles) . The question may then 
be put forward if and why such an extension gives rise to a correct description of irreversible phenomena, as it 
obviously does in the case studied in this paper. The answer is that an extension to internal degrees of freedom is 
legitimate as long as the irreversible processes connected to these variables are slow in comparison to a molecular 
time scale. Think however for a moment that the N particles considered are neither large nor heavy compared to 
bath molecules but form a marked dilute set of these. In that case the Fokker-Planck equation derived makes no 
sense for a description of the motion of the molecules concerned, which is such that the Maxwell velocity distribution 
is attained on a molecular time scale. However, the (slow) long time diffusion regime of the relevant particles may 
still be described by standard nonequilibrium thermodynamics. We have an illustration here of the limitations for 
extending that theory: it is an illusion that a small number of added pseudo-thermodynamic variables, whatever their 
nature , may suffice to adequately describe (with a fundamentally macroscopic discipline) phenomena taking place on 
the time scale of molecular, or microscopic, events. However when one is confronted with a separation of time scales 
and some phenomena take place at intermediate times, a proper extended thermodynamic nonequilibrium theory may 
demonstrate its value in providing the means to easily describe and categorize, together with their symmetries, the 
irreversibilities which take place. 
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